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Abstract 

We estimate the nuclear polarizability correction to atomic transition frequencies in various helium iso- 
topes. This effect is non-negligible for high precision tests of quantum electrodynamics or accurate de- 
termination of the nuclear charge radius from spectroscopic measurements in helium atoms and ions, in 
particular it amounts to 28(3) kHz for 1S-2S transition in 4 He + . 
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There are various corrections which need to be included in an accurate calculation of the atomic 
energy levels. These are relativistic and QED effects, finite nuclear mass and, to some extent, finite 
nuclear size. There is also an additional correction which comes from possible nuclear excitation 
due to the electric field of the surrounding electrons. This effect is usually neglected, as it is 
relatively small compared to the uncertainties in the nuclear charge radii. There are however ex- 
ceptions, where the nuclear polarizability correction can be significant. The first known example 
was the 2P 3 / 2 — 2Si/ 2 transition in muonic helium /i- 4 He + yl Q], where the polarizability cor- 
rection is about 1% of the finite nuclear size effect. Another example is the isotope shift in the 
IS — 2S transition frequency between hydrogen and deuterium. The nuclear polarizability correc- 
tion of about 20 kHz, was two orders of magnitude larger than the experimental precision [3], and 
helped to resolve experimental discrepancies for the deuteron charge radius. Another very recent 
example is the isotope shift in the 2S — 3S transition frequency between n Li and 7 Li Q4£| . The n Li 



nucleus has probably the largest nuclear polarizability among all light nuclei, with a contribution 
to the 2S — 3S isotope shift of about 36 kHz. In this paper we study in detail the nuclear polariz- 
ability correction to atomic transitions for helium isotopes: 3 He, 4 He and 6 He and compare with 
currently available and planned accurate measurements of transition frequencies in helium atoms 
and ions. 

The interaction of the nucleus with an electromagnetic field can be described by the following 
Hamiltonian: 

H int = qA°-d-E-i2-B-l (r 2 ) V • E, (1) 

o 

which is valid as long as the characteristic momentum of the electromagnetic field is smaller than 
the inverse of the nuclear size. Otherwise, one has to use a complete description in terms of 
formfactors and structure functions. Under this assumption, the dominant term for the nuclear 
excitation is the electric dipole interaction. This is the main approximation of our approach, which 
may not always be valid. It was checked however that higher order polarizabilities are quite small 
(below 1 kHz) for deuterium [5, (]], and this should be similar for He isotopes. Within this low 
electromagnetic momentum approximation, the nuclear polarizability correction to the energy is 
given by the following formula |7|] (in units h = c = 1, e 2 = Aix a): 

E P oi = -ma 4 (X^ 3 ( r a)) ( m3 " P oi), (2) 

a 

where m is the electron mass and the expectation value of the Dirac 5 is taken with the electron 
wave function in atomic units. For hydrogenic systems it is equal to Z 3 /tt. In the equation above, 



ctpoi is a weighted electric polarizability of the nucleus, which is given by the following double 
integral 
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where k = a/1 + 2im/w and _&r is the excitation energy for the breakup threshold. The kets | 4>n) 
and \E) denote the ground state of the nucleus and a dipole excited state with excitation energy E, 



respectively. The square of the dipole moment is related to the so called B(E1) function by 

\((f> N \d\E)\ 2 - 



3 dE (4) 
in units e 2 fm 2 MeV -1 , which explains the presence of e 2 in the denominator in Eq. ©. The 
B(E1) function is, in turn, directly related to the photoabsorption cross section at photon energies 
E 

16 7r 3 „ / 1 dB(El)\ 
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which allows us to obtain the B{E1) function from experimental data. 

If Et is much larger than the electron mass m, one can perform a small electron mass expansion 
and obtain a simplified formula, in agreement with that previously derived in [8|]: 
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where a^ is the static electric dipole polarizability and asiog is the logarithmically modified polar- 
izability. We have tested this approximation for 3 He and 4 He isotopes, and found that numerical 
results differ from the exact formula in Eq. © by less than 0.1%. 

In the opposite situation, i.e. when m is much larger that E T , that corresponds to the nonrela- 
tivistic limit, the polarizability correction adopts the form (with m being the reduced mass here): 
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This approximation is justified for muonic helium atoms or ions, because the muon mass (~ 
106 MeV) is much larger than the threshold energy Et [see Eq. ([TO]) !. This formula requires, 
however few significant corrections, namely Coulomb distortion and formfactor corrections. They 
were obtained by Friar in yD for the calculation of the polarizability correction in \i 4 He. This 
nonrelativistic approximation, however, is not valid for electronic atoms since the typical nuclear 
excitation energy in light nuclei is larger than the electron mass. 

We consider in this work three helium isotopes, namely 3 He, 4 He and 6 He, which are stable 
or sufficiently long lived for performing precise atomic measurements. The separation energy 
S = Et for these helium isotopes are Q9Q : 

S 2n ( 6 He) = 0.975 MeV, 
S p ( 4 Re) = 19.8 MeV, 

S n ( 4 He) = 20.6 MeV, (10) 
S p ( 3 Re) = 5.493 49 MeV, 
S n ( 3 He) = 7.718 04 MeV. 

The separation energy S is the main factor which determines the magnitude of the nuclear polar- 
izability correction, since E po \ is approximately proportional to the inverse of S. 

We first consider the 6 He isotope. In this case the polarizability a( 6 He) was calculated from 
two different B(E1) distributions. The first one corresponds to the experimental distribution ex- 
tracted by Aumann et al. [fioll from the Coulomb breakup of 6 He on lead at 240 MeV/u. These 
data are represented by the dots in Fig.[TJ The second B(E1) distribution corresponds to a theoret- 
ical calculation, and was obtained in fl 1 XII by evaluating the matrix element of the dipole operator 
between the ground state and the 1 continuum states. These states where obtained by solution 
of the Schrodinger equation, assuming a three-body model for the 6 He nucleus, with pairwise 
neutron-neutron and neutron-a interactions, plus an effective three-body force. The B(Ei) ob- 
tained in this calculation is represented in Fig. ED by the dashed line. In spite of the discrepancy 
between the theoretical and the experimental B(E1) distributions, the deduced polarizability a po i, 
as obtained from Eq. ©, is similar in both cases: a po i,the = 24.2 fm 3 versus a po i,ex P = 21.1 fm 3 . 
It should be noted that the experimental data were extrapolated and integrated up to Et = 12.3 
MeV, the threshold value for the decay into two tritons. We do not have, however, experimental 
data for the B(E1) beyond this threshold. Therefore, for the final result we take the average and 
add the polarizability of 4 He [calculated in Eq. (TT6l) l to partially account for other decay channels, 




FIG. 1: Electric dipole line strength dB(E\)/dE in units e 2 fm 2 /MeV for 6 He. The solid line is a fit to the 



experimental distribution of Aumann et al. H 1QTI (filled circles), extrapolated to 12 MeV. The shadow region 



indicates the experimental uncertainty of the data. The dashed line is the calculation of Thompson et al. 



[11]. 



and obtain 

a po i( 6 He) = 24.7(5.0) fm 3 = 4.3(9) • 1(T 7 m" 3 , (11) 

where m is the electron mass. For the comparison with other possible determinations of 6 He 
polarizabilities, we additionally present in Table I the static electric dipole and logarithmically 
modified polarizabilities. However, they can not be used to determine a po \ since the Et is of the 
order of the electron mass m. The resulting contribution to the frequency of, for example, the 
2 3 S'i - 3 3 P 2 transition in 6 He is 

E pol = -ma 4 (s 3 ( ri )+5 3 (r 2 )) (m 3 5 pol ) = hv pol . (12) 

is pol = 0.015(3) MHz. (13) 
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Ref. a E [fm 3 ] a E io g [fm 3 ] a po i[fm 3 ] 



3 He 0.153(15) 0.615(62) 3.56(36) 

[12] 0.130(13) 

[13] 0.145 

[14] 0.250(40) 

4 He 0.076(8) 0.365(37) 2.07(20) 

[13] 0.076 

[12 0.0655 

[2] 0.072(4) 

6 He 1.99(40) 4.78(96) 24.7(5.0) 



TABLE I: The electric dipole polarizability a E , the logarithmically modified polarizability aEiog, and the 
weighted polarizability a po \ for helium isotopes. 



For comparison, the finite nuclear size correction to the same transition is I116h 

E is = ^Za 4 m 3 (r c 2 h )(5 3 (r 1 )+5 3 (r 2 )^ __ = hv* 



235i-3 3 P 2 
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-4.253 MHz . 



(14) 
(15) 



The relative magnitude of the nuclear polarizability to the nuclear finite size for 6 He is about 
0.35 %, so it alters the charge radius determination at this precision level. However, the uncertainty 
of the experimental result of Wang et al. ilia! for the isotope shift between 6 He and 4 He, u mo = 
43 194 772(56) kHz, is about four times larger than z/ po i, and therefore the nuclear polarizability 
correction at this precision level is not very significant. 

We proceed now to evaluate the 4 He nuclear polarizability. This has been obtained from the 
total photoabsorption cross section measured by Arkatov et al. DJ/ZL LLE LL9Q- Using Eq. ©, we 
extracted from these data the B(Ei) distribution shown in Fig. [2] (filled circles). Then, applying 
Eq. ©, one obtains the weighted polarizability of 4 He: 



a( 4 He) = 2.07(20) fm 3 = 3.6(4) 10"* m" 



(16) 



It should be noted that at 100 MeV the dipole approximation in Eq. (OQ) may not be correct, since 
the corresponding photon wavelength is of order of the nuclear size. Therefore we introduce 
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FIG. 2: Electric dipole line st rength 



the data of Arkatov et al. I17J, 



18 



dB(El)/dE in units e 2 fm 2 /MeV for 4 He. The circles correspond to 



19|]. The solid line is a fit to the data. 



10% uncertainty to account for this approximate treatment. Results for the static polarizability 
obtained in this work, along with those obtained in other works, are presented in Table I. Our 
static polarizability agrees with that of Friar yfl, which was based on earlier experimental data for 
the photoproduction cross section. It agrees also with theoretical calculations by Leidemann [13], 
but slightly disagrees with the most recent calculations of Gazit et al [ Fsll . 

The obtained weighted polarizability gives a relatively small effect for transitions in neutral 
helium. The ratio of the polarizability to the finite size correction is 4.6 • 10" 4 . At present, there 
are no measurements for the charge radii at this precision level. However, from the planned high 
precision measurement of the 15 — 2S transition in 4 He + |20j], one can in principle obtain the 
charge radius of the a particle from the knowledge of the nuclear polarizability. The polarizability 
correction to this transition amounts to 



u pol (lS - 2S, 4 He + ) = 28(3) kHz 



(17) 



and is smaller than the uncertainty of about 400 kHz in current theoretical predictions, which are 
due to a 2 (Z a) 6 higher order two-loop electron self-energy corrections 112 ill . 
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For the 3 He atom, we separately calculate the nuclear polarizability corrections coming from 
the two- and three-body photodisintegration processes. There are many experimental results for 
the photoabsorption cross section as well as theoretical calculations. They agree fairly well for the 
two-body dissociation, but they significantly differ for the case of three-body disintegration at high 
excitation energies (£7 7 > 20 MeV). Since these experimental results have large uncertainties we 
prefer to rely on theoretical calculations which agree with each other very well. In this work, we 



will use the calculations of Deltuva et al. 112211 and Golak et al. 112311 . The former uses the realistic 
CD Bonn NN potential, supplemented with a three-body force to account for the A excitation, and 
a full treatment of the Coulomb potential. The calculation of Golak et al. is based on the AVI 8 NN 
interaction in combination with the UrbanalX three-nucleon force, and considers Coulomb effects 
just for the ground state. The B(Ei) distributions, extracted from these theoretical photodisinte- 
gration cross sections, are displayed in Figs. 3 and 4, respectively, as a function of the excitation 
energy. 
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FIG. 3: Electric dipole line strength dB(El)/dE in units e 2 fm 2 /MeV for two-body disintegration of 3 He. 



The data are from Fetisov et al. |24l] . Ticcioni et al. 02511 and Kundu et al. 02611 . The solid and dashed lines 

13] and Golak et al. 



are the calculations by Deltuva et al. 112211 and Golak et al. 



respectively. 
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FIG. 4: Electric dipole line strength dB(El)/ dE in units e 2 fm 2 /MeV for three-body disintegration of 3 He. 
The data are from Faul et al. 127D and Berman et al. [28]. The solid and dashed lines are the calculations 



by Deltuva et al. 12211 and Golak et al. 12311 . respectively. 



Using these theoretical B(E1) distributions, the result for the weighted polarizability due to 
both two- and three-body disintegration is 



a( 3 He) = 3.56(36) fm 3 = 6.2(6) l(T 8 m~ 



(18) 



and the related contribution to the frequency of the IS — 2S transition in He + is 



v po i(lS - 2S, 3 He + ) = 48(5) kHz. 



(19) 



Hence, the magnitude of the nuclear polarizability correction for 3 He is almost twice as large as 
for 4 He and can be significant for the absolute charge radius determination from the IS — 2S 
measurement in hydrogen-like helium. The corresponding contribution to the 4 He- 3 He isotope 
shift in the 2 3 £i — 2 3 Pj transition of —1 kHz is at present much smaller than the experimental 
precision of about 30 kHz [29j]. 



Our result for the static po 



arizability of 3 He is presented in Table I. It is in agreement with 



the first calculation by Rinker Ill2ll which was based on measured that time photoabsorption cross 



section, it is also in agreement with calculations of Leidemann [13], but is in disagreement with 
the cross section mesurement for the elastic scattering of 3 He nuclei of 208 Pb below the Coulomb 
barrier Ill4ll . 

In summary, we have obtained the nuclear polarizability correction in helium isotopes. In most 
cases, we find that the correction to the energy levels is smaller than current experimental preci- 
sion, but could affect the determination of the charge radius when more accurate measurements 
become available. Together with possible high precision measurements in muonic atoms, it will 
allow for an improved test of QED and a more accurate determination of the fundamental con- 
stants. 
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